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ABSTRACT 

A unique approach to the problem of gaging propellant tanks 
in a zero gravity environment has been investigated. A capacitance 
technique is used in which surface electrodes are attached to the 
inside tank wall so that the whole tank becomes the inside of a 
capacitor. The shapes and potentials of the electrodes are adjusted 
so that the interior electric field is approximately uniform throughout 
the tank, insuring equal weighting to all regions, regardless of 
how the liquid is distributed. 

A theoretical analysis of the field distribution in a cylindrical 
tank was conducted. This showed that satisfactory uniformity can be 
obtained with a reasonable number of electrodes. It also showed 
that the errors caused by the surface charges on the liquid interfaces 
are not significant in the cases of liquid hydrogen and liquid 
oxygen, particularly if the readings are averaged over two orthogonal 
directions of the field. 

An experimental investigation was conducted in which a cylindrical 
tank with one hemispherical end was constructed, and sixteen ( 16 ) 
electrode strips were installed. A breadboard electronic system was 
built and tested with the model tank. The results were satisfactory, 
except that the method for eliminating the "leakage" capacitance to 
the tank wall was found to be inadequate. New methods for overcoming 
this problem were evolved, but not tested. 
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INTRODUCTION 

With the current interest in manned space missions, a number of new 
technical problems have arisen which must be solved in order to insure the 
safety of these missions. 

One such problem is the accurate gaging of liquid propellants in 
orbiting vehicles, where gravity is essentially absent. Since fluid 
behavior becomes wholly different in zero-gravity, due to the dominance of 
capillary forces, the conventional techniques for propellant gaging are 
obsolete. Hence, new schemes are needed to solve this problem. 

In the gaging of liquid propellants, there are a number of possible approaches. 
Each approach falls into one of two categories, depending on whether the 
propellant mass is inferred from measurements made on the vapor ullage or 
whether the propellant mass is measured directly. 

The ullage techniques have the advantage of being less dependent on the 
geometry of the tank or the orientation of the propellant than are the direct 
mass measuring schemes. The truth of this lies in the fact that they are 
concerned with the measurement of one or more of the thermodynamic parameters 
which are independent of geometry, rather than relying on the measurement of 
some material property (absorption coefficient, dielectric constant, etc.) which 
requires geometry- dependent arrays of sources and detectors. However, the 
ullage schemes are not advantageous insofar as the very thermodynamic parameters 
upon which the measurement depends, in most cases, cannot be adequately 
controlled to be properly measured because of changing thermal inputs, venting 
requirements in some cases, response times, vapor condensation and, in 
general, the persistent presence of thermodynamic transients. Although 
conditions of thermodynamic equilibrium are not necessary in the mass schemes, a 
major drawback is that they are geometry- dependent and require statistical 
averaging techniques to eliminate the resulting errors. 
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Under this contract, a study has been made to determine the feasibility of 
using a mass measuring scheme based on capacitance for solving the zero-g 
gaging problem. The use of capacitance techniques for the measurement of liquid 
level and fluid quantity in pne-g is well established. The proportionality of 
average mass density p to dielectric constant 6 expressed by the Clausius-Mosotti 
relation: 


p = K 


6-1 

6+2 


makes this technique applicable. 


Under the influence of gravity, and in containers of known volume, the ; 
measurement of the effective dielectric constant of liquids of known density 
can be calibrated accurately to read fluid level directly. Similarly, in zero-g 
and in tanks of known volume, a measurement of the effective dielectric 
constant of a two-phase (liquid/vapor) fluid gives a measure of the average 
density and, therefore, the total liquid mass. However, in zero-g, where 
surface tension forces predominate, the configurations of confined two-phase 
fluids are variable and unpredictable. Because these configurations lack the 
usual flat liquid/vapor boundary which occurs under the influence of gravity, 
the conventional point level or continuous level gages are not adequate. Further 
more, the presence of a sensing device within the tank enclosure could affect the 
propellant distribution because of the capillary forces; i.e., probes of the 
usual tubular or coaxial design would tend to draw up and contain the propellant, 
thereby indicating the presence of a greater mass of propellant than actually 
exists , 


In view of these considerations, a completely new capacitance gaging 
concept is necessary. In the determination of a new capacitance design, it is 
necessary to meet spme unusual requirements. The most fundamental of these 
requirements is that the gage must perform in a manner which does not allow any 
region of the tank volume to be weighted more heavily than any other. This 
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requirement arises because in zero-g practically every region (aside from those 
regions near internal structures) is equally likely to contain fluid. Previous 
attempts to meet this requirement have resulted in gage designs comprised of 
a matrix of wire electrodes distributed throughout the tank volume, dividing 
the space into small cells, each monitored independently and then averaged. 
Although this approach seems capable of attaining the desired measurement, it is 
heavy and cumbersome to fabricate. 

Just prior to this contract, however, a capacitance technique was devised 
which circumvents the aforementioned difficulties very well. In the 
present concept, it is not necessary to physically partition the tank space 
into cells in order to monitor each region equivalently. Equal weighting of 
the tank space can be accomplished by establishing a suitable electric field 
configuration within the tank space using electrodes mounted only on the 
tank walls. Furthermore, by properly designing these electrodes, their influence 
on positioning the fluid can be made negligible. 

This new concept attempts to solve the zero-g gaging problem, not 
by incorporating separate probes, but by utilizing the basic tank as a 
"level" gage, whereby the judicious choice of electrodes on the tank walls 
allows the accurate measurement of propellant quantity in the tank almost 
independent of its configuration. 

Background Information : Before discussing the gage concept and the 

results of the feasibility study, it is advantageous, for the purposes of 
reference, to give a brief discussion of the actual SIV-B tank, its internal 
structures, and the physical properties of the propellants and LOX which are 
of relevance in this contract work. 
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The SIV-B Tank : The SIV-B, built by Douglas Aircraft Company, is the 

upper stage of both the Saturn IB and Saturn V rockets. 

The stage is essentially cylindrical and is 58,4 feet long by 21.7 feet 
in diameter. The propellant tank is 21.7 feet in diameter and 44 feet long. 

It features two dome- type ends and an intermediate insulated bulkhead separating 
the liquid hydrogen tank from the oxygen tank. The tanks are made of aluminum 
alloy with the inner surfaces milled in a waffle- like pattern. In the hydrogen 
tank, for purposes of thermal insulation, each recessed "square" is filled 
with a polyurethane and fiber glass block and, when all squares are filled, smooth 
coated with a polyurethane and fiber glass layer. This insulation is about an 
inch thick. Figure 1 illustrates the configuration of the stage and its propellant 
tanks . 

The hydrogen tank contains a number of internal components and structures. 
These include slosh baffles, pressure bottles, instrumentation probes as well 
as the previously-mentioned insulating layer. 

There are two annular slosh baffles in the Ll^ tank. They are made of 
nylon mesh and are anchored to the inner surface of the tank. Each of these 
nylon baffles is about three feet wide and are interwoven with nylon rope for 
support. The first of these is located about ten feet below the tank top and 
the second is further down. There are eight helium pressure bottles attached 
to the inner wall of the tank. They are each spherical with a diameter of 
about two feet. There are also two instrumentation probes which run over 
most of the tank length. These probes hold the sensors for the propellant 
utilization system, temperature measurements, etc. 

The LOX tank has a single centralized instrumentation probe and some 
structural supports and ducting. It does not have the fiber glass insulation 
found in the LI^ tank. 
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Figure 
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Propellant Properties : The propellants considered in this feasibility 

study are the cryogenic fluids LOX and L^. The ratio of density to surface 

tension for these propellants at their boiling point temperatures is of the 

2 

order of a few hundredths of a gm/cm /dyne, which means that even in acceleration 
fields as low as 20 micro-g's, their hydrostatic" behavior will begin to be 
gravity dominated (Bond number^ 1) in tanks as large as the SIV-B. This 
indicates that in low orbit missions, where aerodynamic drag may be of this order, 
the liquid/vapor configuration will approach that in which the vapor is maintained 
"above" the liquid. In a more nearly ideal zero-g environment, the liquid/vapor 
configuration can take a number of forms. However, surface tension, thermal 
gradients, and the almost totally wetting character of these cryogens will tend 
to cause engulfed ullage to migrate along thermal gradients toward the 
warmest portion of the tank wall. Because of the zero contact angle of these 
propellants, the ullage will form nearly spherical interfaces with the liquid. 

However, these simple and completely static configurations will be 
constantly perturbed by external accelerations caused by attitude control, 
venting, etc., and by propellant boiling and cavitation. For these reasons, 
the chosen gaging scheme and the experimental portion of the feasibility study 
does not rely heavily on these configurations but attempts to evaluate the system 
under more general considerations. 

The electrical properties of LH^ and LOX are highly appropriate for the 
present gaging technique. These propellants have extremely high electrical 
resistance So that their electrical response is completely reactive, which is 
required. Also the dielectric constants are small enough so that errors introduced 
by dielectric shielding are not large, yet these constants are big enough to 
allow a substantial change in gage impedance between empty and full tank. The 
relative dielectric constant of is about 1.23 and that of LOX approximately 

1.5. 
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PROGRAM PURPOSE 

This report describes the work that has been performed under NASA 
Contract NAS8- 18059. The total contract period covered is from July 25, 1966, 
to January 15, 1967. 

Contract Objective : The objective of the contract has been to perform 

a feasibility study of capacitance techniques for mass gaging of propellants 
in zero-g. In particular, the chosen technique was to be capable of performing 
to all requirements listed below when used in the main LOX and LH^ tanks of 
the Saturn SIV-B stage, as well as in small tanks on the order of 8' x 4'. 
Emphasis was to be placed on proving feasibility for the SIV-B stage, since it 
was assumed that this was enough to insure feasibility for smaller tanks. 

The basic performance requirements of the gaging system were set forth 
as follows: 

The propellant measuring system shall measure the total mass 
of propellant (liquid hydrogen or LOX) in the tanks to an 
accuracy of ±2% of full scale or better. 

This accuracy shall be achieved only when the propellant is 
essentially in static equilibrium such that the vapor ullage 
is confined to a singljb ’’bubble" located somewhere in the tank. 

This accuracy shall be achieved in acceleration environments of 
zero-g to five g's. However, if special considerations are 
necessary to achieve this accuracy in higher g's, emphasis shall 
be placed on meeting this requirement in zero-g. 

The response time of the propellant gaging system shall be 0,5 seconds 
or less. 

To meet these requirements, a new capacitance scheme was devised and 
evaluated under this contract. 
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The tank was to be designed in all essential details including 
mock-ups of internal structures (slosh baffles, helium bottles, probes, etc,). 
Furthermore, scaling laws for the model were to be established so that a 
valid prediction of the performance of a full-scale gage could be made from 
the test data for the model. 

For the purposes of evaluating the experimental performance of the model 
gage, a zero-gravity environment was to be simulated by using either two 
equidensity, immiscible fluids with different dielectric constants or by 
simulating zero-g fluid configurations by appropriately molded materials of 
low dielectric constants (such as polyurethane foam). 
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I. GAGING SYSTEM CONCEPT 

The capacitance concept studied under this contract is designed to meet 
certain physical requirements which are necessary to insure the development 
of an accurate zero-g propellant gaging system. The underlying effort 
in the approach is to devise a design which takes full advantage of the 
tank geometry in the generation of the sensing field so that unnecessary 
electrode structures within the tank space are avoided. 

For this reason, the concept utilizes the tank walls as a desirable 
location for the placement of the source and sensing electrodes. The 
specific physical design of these electrodes and the potentials applied 
to them are determined by both the tank shape and the desired field 
configuration within the tank space. 

1.1 Ideal Field Configuration : In a zero-g application, it is necessary 

to generate a field which samples each point of the tank space 
equivalently so that liquid in one region is not weighted more 
heavily by the measurement than an equal amount of liquid elsewhere. 
Since*the weighting is proportional to the field strength at each 
point, this is accomplished best by a field configuration which is 
uniform; i.e., identical to the field of a parallel plate capacitor. 

To generate a uniform field in a tank which has curved boundaries, 
it is required that a distribution of source potentials be impressed 
on these boundaries. In fact, to generate a perfectly uniform field 
inside of a tank with curved walls, this potential distribution must 
be a continuous function of the coordinates specifying the wall 
geometry. To be more specific, in Appendix A we have calculated the 
required boundary potential functions for a cylinder and a sphere 
when the uniform field is es tablished transverse to the tank axis,* 

*It is also possible to choose the direction of the field parallel 
to the tank axis; however, this is unattractive for the following 
reasons: First, the electric field must be parallel to the cylindrical 

portion pf the tank and this would require a considerable number of 
electrodes (circular bands) to maintain. Second, the boundary potential 
must be a linear function along the tank length. This would mean that 
none of the individual potentials would be the same, as in the present 
case, so that considerably more transformer taps would be necessary. 
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For example, at the cylindrical wall of the LH^ tank, the 
circumferential variation of potential must follow the relation: 

V (0) = V q cos0 (1.1) 

while for both the hemispherical domes of the LH^ tank and the entire 
LOX tank (assumed spherical) the potential must obey: 

V (0, 0) = V q sin 0 cos0 (1.2) 

The value is determined by the product of the desired electric 
field E and the tank radius R, and these expressions are the correct 
forms which match E and V along the dome-cylinder junctions. 

As a consequence of the required variations of potential along the 
cylindrical or spherical tank boundaries, it is not possible to have 
a single conductor, like the tank wall, possess this boundary potential. 
For a real gage, the ideal case can be approximated only by incorporating 
within the tank an array of separate conductors which conform to the 
geometry of the walls and which carry potentials in accordance with 
the relations (1.1) or (1.2), depending on the angular coordinate (s) 
of the conductor position. For example, on the cylindrical wall, a 
number of conductor strips of prescribed width extending the length 
of the wall can be bonded to the wall insulation (or fastened close to 
the wall by other means) so that each strip is parallel to its neighbor 
and separated from it by a small distance. To each strip the correct 
potential is applied using (1.1), Obviously, the more strips that are 
included, the more exact is the approximation to the ideal case. 

However, if a very large number of strips were required, it would not 
be a very attractive scheme for mass gaging. It turns out, as we have 
theoretically and experimentally verified and discuss in Section III, 
that a large number of electrodes is not required and that fewer than 
twenty electrodes produces a field of considerable uniformity. 
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Applying the game technique to the spherical LOX tank or the 
hemispherical end caps of the tank, the problem is more complex 
because of the required potential variation in the two angular 
directions, Q and 0. Here, use can be made of a number of 
"patch sensors" which are excited with potentials in accordance with 
Equation (1.2). 

Because of the necessity of using a discrete array of boundary 
electrodes in a real gage, it is not possible to generate a 
perfectly uniform field. For this reason, one of the important 
tasks in the feasibility study is the determination of the number, 
widths, spacings, and angular positions of the conductors in 
the array to produce a field of optimum uniformity for the application. 
We present the electrode design analysis in Section III. 

1.2 Averaging Concept : There is a second physical requirement which the 

gage concept may have to be capable of meeting. This requirement has 
to do with the dependence of the gage capacitance on propellant 
orientation. Even though the use of a uniform field guarantees that 
equal fluid elements at any point in the tank space will contribute 
equivalently to the gage capacitance, strictly speaking, this is only 
true if the fluid elements have the same shape and orientation. This 
effect arises as the result of the induced surface charges which the 
field creates at dielectric interfaces. The magnitude of these 
charges depends on the direction of the field relative to the 
boundary. Fluid elements with surfaces parallel to the flux lines add 
directly to the dielectric value (parallel capacitance case) while 
elements with boundaries normal to the field lines always contribute 
less (series case). The effect is known as dielectric self- shielding 
and is greatest for materials of high dielectric constant. In the 
present application, the fluids LH 2 and LOX have fairly low dielectric 
constants so that errors introduced by shielding will be correspondingly 
small. However, it is possible to reduce these errors by an averaging 
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technique which involves applying two orthogonal fields on a 
time sharing basis, or by completely rotating the field about the tank 
axis and averaging the individual measurements. It is not clear at 
this point whether averaging is really necessary with LH 2 or LOX 
since their dielectric constants are small. Furthermore, if the 
supposed spherical ullage bubble configuration is actually valid in 
zero-g, the averaging scheme would not be effective anyway, because 
of the bubbles' symmetry. For these reasons, and because the 
inclusion of this technique would complicate the electronics 
considerably, the experimental portion of the present program did 
not include averaging. Instead, stress was placed on evaluating 
the basic physical and electronics concepts to determine if a 
simpler design could achieve the desired characteristics. The need 
for averaging would then depend on whether or not the simpler 
concept could be made adequate. 

1.3 Measurement Concept : The final requirement which the gage concept 

must meet is that the measurement technique must result in an output 
signal which is proportional only to the quantity of propellant 
in the tank. Furthermore, this signal must be large enough so that 
the minimum signal resolvable by the measuring circuit is well below 
it and within the desired accuracies. 

The measurement problem is illustrated by the simplified circuit 
shown in Figure 1-1. The problem is to measure the capacitance C, 
which is a function of the propellant quantity, while maintaining 
electrodes 1 and 2 at potentials dV and -V respectively. These 
potentials are necessary in order to provide the approximately uniform 
field so that the measurement will be relatively independent of the 
propellant distribution in the tank as described in an earlier 
section of this proposal. 
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The measurement must be made by a scheme in which the result is 
not affected by variations in Cj and C 2 which are large compared 
to C. C-^ and C 2 consist of capacitances from the electrodes to 
the tank wall and cable capacitances. C-^ and C 2 are of the order of 
1000 times the full-scale value of C. 

Finally, the measurement scheme should "subtract out" the 
empty- tank value of C. The capacitance C may be written: 

csc o [i + (e-D 5 ] d.3) 

where C q is the value of C when the tank is empty 

£ is the relative dielectric constant of the liquid propellant 
m is the mass of liquid propellant 

M is the capacity of the tank in terms of liquid propellant mass. 

The quantity to be measured is the second term; i.e., 

Ac = c - c K (e - l) § c (1.4) 

o Mo 

which is a measure of the propellant mass. 

The circuit of Figure 1-2 performs the desired measurement. The 
operation of this circuit may be described more simply by the use of 
an equivalent circuit referred to the primary side of T^. This circuit 
is presented in Figure 1-3. The system is excited by an oscillator 
which provides the potentials V and -V, which are alternating at a 
low RF rate, for example, 100 KHz. With the switches in the 
position shown, 

I = V ( jwC, + jwC + jw4C~ jw4C ' +— — ) 

1 1 2 0 jwL x 
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Figure 1-2 
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and with the switches in the alternate position. 


I 2 - VOwCj + jwC 2 +'- 55 -) 


The difference between these is. 


I x - I 2 = V ! jw (4C - 4C' o )] 


(1 = 5) 


If C' is adjusted to be equal to C , the empty-tank value of C, 
then 1^ - Ig is proportional to the AC of equation (1.4). 

The current I is, therefore, a square- wave-modulated signal in which 
the peak-to-peak value of the modulation is proportional to the 
propellant quantity. Although the "carrier" current is a function of 
Cp C 2 > and L^, which are large quantities and may not be stable, 
the modulation is independent of all of these quantities. 


In order to have an accurate measurement it is necessary that V 
and CO , the amplitude and frequency of the source, be stable to somewhat 
better than the measurement accuracy, and also that the amplifiers and 
detectors operating on E be similarly stable. Stabilities on the 
order of 1/4-7° are required; such stabilities are readily achieved in 
reliable space- tested hardware. 


In the actual measuring system, the transformer T^ has a secondary 
with multiple taps in order to provide the proper potentials to the 
various electrodes. Figure 1-4 shows the system for a twelve-electrode 
configuration, for one of two spatially orthogonal fields. With the 
transformer winding as shown, the capacitance C q is approximately 


C 

o 
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where 6 is the absolute dielectric constant in farads per meter 
and X is the length of the tank in meters. For a tank of 10 to 15 
feet long, C q is approximately 30 picofarads. Of course, by 
altering the turns-ratio of the transformer, this capacitance 
(referred to the primary) may be made any convenient value. To make 
the system less sensitive to propellant distribution, two spatially 
orthogonal fields are used on a time-shared basis, and the outputs 
aye averaged. The two systems utilize the same set of electrodes 
within the tank. Semiconductor switches operated at a low frequency 
(e.g., 500 Hz) alternately connect the electrodes to the potentials 
as shown in Figure 1-4 and in a similar configuration rotated 90° in 
space. 
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Capacitance Mass Gaging System 
Figure 1-4 
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II. PERFORMANCE ANALYSIS OF AN IDEAL GAGE 

As a preliminary to the theoretical analysis of the actual gage, we 
present here a performance analysis of an ideal gage in order to 
demonstrate the basic soundness of the uniform field concept. The approach 
here is completely analogous to that used in a later section where the 
analysis of the actual gage is presented. However, the mathematical 
treatment in this case is much simpler and as a result, more lucid. 

We will perform two basic calculations. The first will be a calculation 
of the approximate empty tank capacitance of the gage to show that this 
concept is capable of meeting the requirements of the measuring circuit; 
i.e., the minimum capacitance criterion. The second calculation will treat 
gage performance when propellant is included in the tank to determine 
whether or not the concept can perform with the desired ±2% accuracy. 

The assumptions which are made in this analysis are the following: 

First, the gage is assumed to be ideal. This simply means that it represents 
the ideal limit of the actual gage so that a continuous cosG potential 
can be established on the tank boundary. This implies a gage with 
infinitely many boundary electrodes, each electrode strip of infinitesimal 
width. Second, the hemispherical endcaps of the tank are neglected and 
the tank is viewed as a cylinder of infinite length. This allows a simpler 
solution of LaPlace's equation. The capacitance calculated from this 
model is, therefore, one per unit of tank length. Third, in the performance 
calculation, the error is computed as a percent (of full tank) deviation 
from the ideal linear response. 

2.1 Tank Capacitance : Aside from the errors introduced by stray 

capacitances, propellant orientation, etc., there will also be some 
intrinsic error associated with the measuring circuit. This error will 
place a limit on the resolving power of the circuit for measuring 
the gage impedance. This, in turn, will require that the gage possess 
some minimum (empty tank) capacity in order that this intrinsic error 
is well within the desired gage accuracy. 
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If we designate the empty- tank capacitance of the gage as C q , 
then the change in gage capacitance from empty to full tank is 
given by. 


AG = ( €/ e -1) C (2.1) 

o o 

where 6/ q q is the relative dielectric constant of the propellant. 
Now, if the limit of resolution of the measuring circuit is &C, 
then this will correspond to some intrinsic volume error AV in 
measuring propellant quantity, expressed by: 


AV - Zc V o 


( 2 . 2 ) 


where V is the total tank volume. Therefore, since it is desired 
o 

to keep AV small compared to errors contributed from other sources, 
a condition is implied on C q ; i.e., 




(2.3) 


Presently, Trans- Sonics is producing quantity gaging circuits 

-3 

capable of resolving capacitance changes of the order of 7 x 10 
picofarads. Using this value, assuming LH^ is the propellant 
(worgt case = 1.23) and requiring that the intrinsic quantity 

error AV/V q should not be greater than 0.1%, the minimum empty 
tank capacitance must be 

C « 30 pf 
o r 

The present gage concept is capable of meeting this requirement. In 
the next section, the expression for the cylindrical tank capacitance 
is derived and, rewritten here, is 

C o = f G o L ’ (2 * 4 > 
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- 12 

where L is the tank length in meters and “= 8.85 x 10 f/m 
(vacuum value). This gives a total empty tank capacitance in the 
vicinity of 50 picofarads or more for the S-IVB tank. 

Since the change in capacitance from full tank to empty tank for 
LH2 will be about 207> of this, or lOpf, it is expected that mass 
changes of the order of a 0.17 o of full scale will be resolvable. 

Thus, it is clear that no difficulties arise with respect to 
signal to noise when resolving mass quantity to within one or 
two percent. 

2.2 Performance Analysis : A derivation of the error signal as a 

function of dielectric constant and propellant quantity is given 
below and is calculated for the case in which the propellant/ 
vapor interface is a cylindrical surface concentric with the tank 
wall. The fuel is contained in the region between this interface 
and the wall. See Figure (II- 1) . This configuration can be 
considered as a gross idealization of the zero-g situation in which 
. the propellant is capillary dominated and completely wets the 
wall, engulfing its ullage. In a real situation, the ullage bubble 
would migrate to the wall, usually to a warmest point. However, 
for our purposes we can obtain from this case an idea of the error 
due to fluid configuration. Furthermore, it is an interesting case 
in that its circular symmetry does not make it amenable to any type 
of field averaging (as would be the case with a spherical bubble 
in an actual tank); and, therefore, the errors are irreducible 
except by calibration. 

Consider a cylindrical tank of radius R with its center at the origin 
of a cylindrical coordinate system. Assume the potential 

V (R) = V q cos© (2.5) 
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as the given boundary condition. Now assume that the boundary 
between a liquid of dielectric constant on th e tank walls and 
a gas of dielectric constant 6^ in the center of the tank is 
a cylinder of radius a, coaxial with the tank. 

The solution of Laplace's equation gives the expression for the 
potential in the gas as 

= Arcos© (2.6) 

and that in the liquid as 

= (Br + c/r) cos© (2.7) 

where A, B, and C are constants to be determined by considering 
the field conditions at the dielectric boundary at r = a. These 
conditions are 


and 




V x (a) = V 2 (a) 


( 2 . 8 ) 


Tha t is. 


e^cos© = £ 2 (B-C/a 2 )cos© 


(2.9) 


and 


A a- = B3- + c/a 


( 2 . 10 ) 
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At r = R we have, from (2.5) and (2.7), 

BR + C/R = V q (2.11) 

Solving (2.9), (2.10), and (2.11) yields 


lo 

1+P 

R 

1+P 2 p 


1 

R 

l+P 2 jS 

Vo 

2 

R 

1+ P 2 jS 


where 

a ill 

and 

p = a/R 


A capacitance per unit length may be defined as the charge on 

one half of the cylinder, per meter of length, divided by the 

potential 2V ; thus 

° tt 


C = 



” z 

Re l (e>-c/R z )cosed& 

J.TT 


or 


7. 



I ~ P 2 ft 

- 1 + P> - 


( 2 . 12 ) 
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If C q is the capacitance of the empty tank; i.e., for a = R, then 

C = -7 6 - and 

o 4 1 


C 




“ 6 1 >- 


_ 2 2 ^ 

1+P 2 |S _ 


(2.13) 


Now, if (1 - p^) = v = normalized liquid volume, (that is, v = liquid 
volume divided by total tank volume), then we would like to have 


AC 


ideal 



( e 2 " * 


(2.14) 


AC - AC 
AC 


max 


ideal 


v (1-v) 
1 +Q-v) 


(2.15) 


The results of this calculation are summarized in Figures II-2 
and II- 3. In both figures, £>C is the capacitance deviation from 
linearity and AC is the total change in gage capacitance from 
full to empty tank. In Figure II-2 the fractional Sc/AC error is 
given as a function of the fraction of propellant in the tank. 

Figure II- 3 displays the maximum error (maximum points of Figure 
II- 2 ) as a function of relative dielectric constant of the fuel. 

Just for reasons of comparison, the error bands for several 
fuels have been computed. It is interesting to note from (2.15) 
that even in the case of 6 -»oo(See Figure II- 3) the worst deviation 
is about 17%. For the fuels of interest, hydrogen and oxygen, the 
maximum uncalibrated errors are about 2 % and 4% respectively and 
this occurs only when the tank is about half full. Calibrated, these 
errors become ±1% and ±2% respectively. Since these errors are based 
on an unlikely fluid configuration as well as a condition which 
does not involve averaging, it appears clear that the gage concept 
proves quite feasible for measuring propellant mass to within the 
desired accuracies. 
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III. GAGING SYSTEM DESIGN 

The design of the actual gaging system is concerned primarily with 
two areas. These include the design of the electrode configuration 
and the design of the measuring circuit. The ultimate design of the 
measuring circuit depends, for the major part, on the resultant electrode 
configuration. For this reason, we will consider the circuit design 
last. 

The design of the electrode configuration is of utmost importance in 
determining the ultimate achievable accuracy of the scheme. For 
example, in the previous section we demonstrated that the desired accuracies 
could be achieved if a perfectly uniform field is generated within the 
tank. However, even with the ideal gage configuration, the achievable 
accuracies (at least for the propellant configuration analyzed) do not 
leave a large margin for additional errors. For this reason, the 
ability to generate a highly uniform field by a finite set of discrete 
electrodes on the tank boundary is crucial to the success of the concept. 
Thus, in order to ascertain concept feasibility, the electrode design 
study is necessarily comprehensive; and is, therefore, a major task of 
the program. 

The analysis performed during this contract, in determining the optimum 
electrode design, is extremely involved mathematically and contains a 
considerable volume of tedious calculations, both formal and numerical. 

To avoid cluttering the main issues of the analysis, the primary attempt 
in this section will be to present the fundamental theory underlying 
the mathematical approach. As will become evident, much of the initial 
computational effort is now less important since the conclusions drawn 
from it are sufficiently simple, in retrospect, that some of these 
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calculations seem superfluous. However, at the time, this was not 
so because the conclusions drawn could not be known a priori, and it 
was only after these calculations were performed that the picture became 
adequately clear. 

Thus, for the purposes of this report, only the main conclusions will 
be presented here. However, sufficient formulation will be included in 
order that the discussion of this analysis is better self-contained 
and of value, if a comprehensive review of the analysis is desired. 

For details of particular calculations, reference can be made to either 
the appendices of this report or the progress reports, 

3.1 Electrode Design Analysis : In designing the electrode configuration, 

there are two considerations to keep in mind: 

1. The electrode configuration must be capable of generating 
a highly uniform field. The degree of uniformity 
required will depend on the gaging accuracy desired, and we 
saw as a result of the analysis of the ideal case that the 
degree of uniformity must be high. 

2. The electrode configuration should conform to one of a class 
of symmetrical arrangements on the tank boundary, which will 
allow the use of field rotation techniques to reduce the 
effects of dielectric shielding, if required. 

These two considerations specify the type of problem to be solved. 
The first consideration requires that the physical characteristics 
of the electrodes (number of electrodes, their sizes, shapes, and 
orientations) and the potentials applied to them, must be chosen 
in a way that results in optimum field uniformity. 
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The second consideration limits this choice by imposing certain 
restrictions on the possible electrode shapes and orientations. 

For example, if complete field rotation about the tank axis is to be 
used, then each electrode must be of the same width, the width being 
dependent on the number of electrodes used. The number of electrodes 
must be even and they must be distributed about the tank circumference 
with equal spacing. If simple averaging in only two (orthogonal) 
directions is to be used, then the electrodes must be disposed symmetrically 
about the 45° axis between the two field directions. For this case, only 
certain pairs of electrodes must have equal widths. However, by an 
appropriate arrangement of the electrodes in the first case, the 
symmetry of the second case can be achieved as well. Therefore, in the 
design analysis, the condition that all electrodes have equal widths is 
imposed. 

It should be mentioned at this point that if averaging is not to be 
used, then the equal width condition is an undesirable one because, for 
a given number of electrodes, it does not lead to a field of optimum 
uniformity. For a field which is to be held fixed in one direction, the 
best uniformity is obtained when the electrodes have decreasing width 
with decreasing magnitude of potential, so that the electrode array 
has 180° symmetry about the tank axis. 

In view of the fact that a fixed field, however uniform, cannot compensate 
errors introduced by shielding effects, it was deemed wise to design the 
electrode configuration assuming that field averaging might eventually have 
to be used. This was done in spite of the restriction it imposes on the 
initial field uniformity. In essence, a practical approach was taken in 
which the resulting design would be an optimum trade-off among field uniformity, 
insensitivity to propellant orientation, and simplicity of practical 
hardware , 
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3.2 Mathematical Model : The major step in the design analysis was 

the determination of a suitable mathematical model of the gaging 
system which was conducive to theoretical analysis and which would 
ultimately lead to a definition of the basic electrode configuration. 

At the outset, it was hoped that the actual system could be 
analyzed in detail without approximations; however, a preliminary 
analysis ^ showed that this was analytically difficult and, therefore, 
an unattractive approach. 

The realization of this fact led us to develop an approximate 

two dimensional model which turned out to be amenable to 

conventional techniques of analysis. In the establishment of 

a tractable model, it was borne in mind that the basic features 

must approximate the actual case very closely in order to determine 

the design parameters accurately. Furthermore, an accurate 

model was necessary to insure a valid performance evaluation of the gage 

when used in analyzing the effects of propellant distribution, 

shape, etc., later in the program. The approximations made in 

establishing the model turned out to be extremely valid; and, as 

a result, the model is an excellent representation of the actual 

case . 

As a result of this fortuitous situation, a calculus of variations 
technique could be used to determine the required electrode sizes, 
positions, and potentials, which are necessary to produce internal 
fields of optimum uniformity. In fact, the general case of N 
electrodes was completely solvable. This was the outcome hoped for, 
and it led to closed form expressions for the desired design parameters 
(electrode sizes and potentials). 


(1) Monthly Progress Report; August, 1966, 
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Because of the basically cylindrical nature of the tank, the two 

dimensional model depicted in Figure (III-l) was determined to be 

the desired one. The model differs from the actual tank in two 

ways. First, this particular model (tank cross-section) implies 

that the tank is a very long cylinder and neglects the end caps. 

Second, the model does not include the insulating layer between 

( 2 ) 

the electrodes and the grounded tank wall. It can be shown 
that the insulation does not appreciably influence the field inside 
the electrode boundary since its thickness is small compared to the 
tank radius. In three dimensions, the electrodes are long, 
parallel conducting strips of equal widths. 

3.3 Optimization of Field Uniformity ; The method used for determining 
the desired discrete electrode configuration which produces a 
field of optimum uniformity is based on a standard calculus of 
variations technique. The usefulness of this technique is made 
clear by the following considerations: 

Let us assume that we place N equal width (angular width = 2 (X) 
equally spaced electrodes on the tank boundary (a configuration 
which would allow averaging) and apply to each of these electrodes 
a fixed potential V , which may be different for each electrode. 

As an example, the distribution of electrode potentials at the tank 
wall might look like that depicted in Figure (III-2). Now this 
potential will establish a unique field pattern within the tank. If 
we change any of the potentials V or change the widths or spacings 
of the electrodes, the field pattern inside the tank will experience 
a corresponding change. The object of the design analysis, therefore, 
is to vary (mathematically) all of these electrode parameters until 
the field pattern becomes the most uniform. For a fixed number of 
electrodes on the boundary there will correspond a best electrode 

(2) Loc. Cit, 
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width, spacing, and potential distribution which results in 
optimum field uniformity. Clearly, this uniformity will become 
more nearly perfect as the number N of electrodes is increased, 
since in the limit of an infinite number of electrodes (with a 
cos 0 potential distribution) the field becomes perfectly 
uniform. 


The mathematical technique for handling an optimization problem 
of this type is the following: First, we analytically determine 

the potential function $ N (r,0) that would exist inside the tank 
due to N equal width electrodes symmetrically positioned on the 
tank wall. This determination requires the solution of Laplace's 
equation (polar coordinates): 


V 2 0. 


N 


9r K dr ’ r 


= 0 


(3.1) 


Inside the circular region of Figure (III-l), subject to the 
N-electrode boundary potential profile V(0) illustrated in 
Figure (III- 2). The general form of (r, 0) for these cases 
will be 


0N ( r > 0 ) 


co 


= ! (a^ r m sin m0 + B m r m cosm©). 


msl 


(3.2) 


where the coefficients A m and Bj^ are determined by the 
particular boundary potential function V(0) through the usual 
Fourier relations, 



^ rram 

r 

l V(0) sin m0d0 

(3.3) 

and 

B L. 

7f a m 

s 

J-TT 

r * 

1 V (0) cos m0d0 

(3.4) 


-'-TT 



where a is the tank radius. 
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Now V(0) contains all of the information about the electrode 

configuration; i.e., V(0) is also a function of electrode width 

206 , the electrode positions and the individual electrode 

potentials Therefore, the coefficients Ajjj and B m also contain 

this information through (3.3) and (3.4). The optimization 

problem, then, becomes one of adjusting the coefficients and 

B in Equation (3.2) by mathematically varying each V and the 
m. is. 

widthsodso that 0^(r, 0) becomes, in functional form, as close 

to the perfectly uniform field potential function 0 q = V q — cos0, 

as is possible. This amounts to saying that we want to adjust the 

coefficients A and B so that the difference A0 = 0 (r, 0) - 
mm o 

0 N (r, 0) is as small as possible at each interior point (r, 0) of 
the tank. 

To do this, we utilize the Mean Square Error Integral Function (3) 
defined by 

r a 2 

M(°6, Vr) = ^p~2 £ 0 o (r,0)-0 N (r,0, <* ,V fc ) rdrd0. (3.5) 

J-TT J 0 

This integral gives the squared difference of 0 q and 0^ at each 
interior point (r, 0) averaged over the domain of interest. 

By minimizing this with respect to the variable design parameters 
and in 0^; i.e., 

i|— = 0; = 0 (for each k) (3.6) 

k 

we can determine the values of and which give the greatest 
uniformity. In fact, the method is completely analogous to a 
"least squares fit" of 0^ to 0 q , The calculation of this integral 
for specific electrode configurations is the primary concern in the 
electrode design analysis, 

(3) See Appendix B for a discussion of this technique. 
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To calculate the Mean Square Error Function of Equation (3.5) we 
need the explicit form of the interior potential 0 (ot,V , r, 9), 

3.4 Determination of 0^: In this section, we present the results of the 

derivation of the interior potential function 0^(r due to 
N electrodes symmetrically positioned on the boundary r = a , Each 
electrode is assumed to have an angular half width, OC . The electrodes 
occur in pairs; one with applied voltage H-V^j the other with -V^. As 
we will see, the magnitudes of the voltages approximate a cosO 
variation around the boundary. 

Before writing down the results for 0 , a point concerning the 
symmetry of the electrode configuration should be made. Since the 
electrodes must occur in pairs, the number of electrodes ;(N) must be 
an even integer. Now these can be broken down into two distinct 
symmetry groups; one designated, 

N = 4k + 2; k = 0, 1, 2, 3 ..... 

“ 2, 6, 10, 14 ..... 

and the other, 

N=4k; k = 1, 2, 3 ..... 

=4, 8, 12, 16 ..... 

The first group has field rotational symmetry, but not 45° symmetry. 
The second possesses both. In addition, the second group can be 
further divided into two subgroups. The distinction here is one 
of electrode orientation; whether a gap or electrode appears at the 
boundary point 9 = The subgroups can be transformed into each 
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other by rotating the electrode system through an angle “. In the 
derivations of 0^ both groups were considered, resulting in an 
expression 0^ for each case. 


The appropriate N electrode boundary potential profile for the N = 4K+2 
symmetry group is depicted in Figure (III-2), and the one for the 
N = 4K group is shown in Figure (III- 3), 


The solution of Laplace’s equation (3.1) inside the region of 
Figure (III-l), subject to these boundary conditions, leads to the 


following explicit forms for 0 (r, 0,C£, V ), For N = 4K + 2: 

N K 

60 • isj 


$/w,v k )= ( v ° + cos ~r) 

tvt-1 k .| / 


cos m© 


(3.7) 


For N = 4K: 


V / V 


^fee,«,V K )= [v„+z £ V.cos 


Sin moL 
m 


COSm© 


(3.8) 


In both of these expressions, the index m takes on only odd values. 


It is to be noted that the only difference in these two expressions 
is in the upper limit of the k sum. Therefore, in the remaining 
calculations we will not distinguish between the two cases. The 
only time that this is necessary is when an explicit writing out of 
the expressions is required. We will not have occasion to do this 
in the following discussions. 

3.5 Calculation of the Mean Square Error Integral Function : With the 

expression for 0^ given by (3.7) or (3.8), and using the uniform 
field potential 0 (r, 0) = V a — cos 0, the Mean Square Error Integral 

Function (3.5) can be calculated. Without doing the tedious computation 
in detail, the result is: 


- 29 - 



TRANS-SONICS, INC., Burlington, Massachusetts 




m <‘Vk 


) - V o 2 j - | AjCKjsinK + j y 


. ... ± / .2 sin mot 

(N) sinOC + — / A m OO^^j 


(3.9) 


rnsi 


where the coefficients A^ (N) and (N) are given by the 
complicated series expressions 

k ( 12 


A. (N) = 

l rr 


v a 


+ 2 


Z 




Vk 2KTT 

V a C ° S N 


(3.10) 


and (introducing new summation indices to avoid confusion) 


A Z (H) s — 





E 

r*i 


Vr 


cos 


2mrTf 

N 




<LL 


V s Vt 

Vi 


cos 


2m sir 
N 


cos 


■zmtrr 

N 


(3.11) 


In all of the sums appearing in these coefficients, the upper 
limits are either or as i- n the 0^ expressions. Because of 
the complicated expression for M (e>6, v ) we will not attempt any 
explicit evaluations. In the actual analysis, each of the trigono- 
metric series appearing in M were evaluated in closed form, so that 
numerical evaluation of the series could be performed easily. However, 
these evaluations are not necessary here as we will see below. 


3.6 Variational Procedure for Determi ning oi and Vjr: To determine the 

optimum design parameters oC and V , the expression for M appearing 
in (3.9) must be minimized, as we have stated; i.e,, we must compute: 
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This leads to a set of simultaneous equations in the unknowns 

and V , the solutions of which give the desired values of these 
parameters. Without details, we write out the resulting simultaneous 
equations : 


For f~ = 0, 




2t T)Ct 


''Wrt 


0 


= o 


„ 3M A , r; UTTAIO N “ 2 N 

For — = 0, and for each K = 0, 1, 2 —7— or -7. 

O V|^ 4 4 

QO 

cos Jgt sm «. _ 2 y C05 2 ^ . , o 


W1S| 


(3.12) 


(3.13) 


The set of simultaneous equations (3.12) and (3.13) are the 

important relations for determining the required design parameters 

and V . The error function (3.9) is the expression by which we are 
K 

able to judge how uniform the resulting field will be; and, therefore, 
it is used to determine the number of electrodes necessary to 
achieve a desired uniformity. 

3.7 Design Analysis Conclusions : If equations (3.12) and (3.13) are 

solved for the cases of N = 2, 4, 6, etc., an unphysical solution 
occurs for the electrode width (X in all cases of N greater than 2. 

The achievement of a minimum in the mean square error requires, in 
these cases, that the electrodes be slightly overlapped. However, 
since this is physically impossible, the conclusion to be drawn from 
this is that M cannot be completely minimized. The smallest value 
of the error function is obtained when all electrodes are butted 
edge to edge; i.e.. 


0 C = tT/n 


(3.14). 
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Furthermore, the corresponding condition on the applied potentials 

V can be shown to be: 

K 


V Tr = V cos 
K o 


2KT r 

N 


(3.15) 


Therefore, the first conclusion of the theoretical design analysis 
is that, for a given number of electrodes, optimum uniformity 
is achieved if the entire circumference of the tank wall is covered 
by the electrode array, leaving very small interelectrode spacing. 
The second conclusion, which is a bit surprising, especially for a 
small number of electrodes, is that the potentials which are to be 
applied follow the V^cosQ relation exactly; the proper value of V^. 
being determined by the angular position of the center of the 
conductor strip. 


Although these conclusions seem evident in retrospect, they could not 
be obtained justifiably with any degree of certainty before the 
calculations were made. For example, instead of computing the mean 
square error over the volume of the tank, one could have computed it 
over the tank circumference, in an attempt to fit a discrete step 
potential function most closely to the continuous cosQ one at the 
boundary. A calculation of this leads to the same conclusions 
concerning the butting of electrodes but the conditions on the V K 
are significantly different. An analysis of this difference 
leads to the conclusion that the second choice results in an optimization 
of field uniformity near the tank wall but results in a bad field 
uniformity in the tank center. In a following section, we will show 
that the volume error approach is the correct one since the field 
uniformity throughout the tank increases at every point with an increase 
in the number of electrodes, as it should. 
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3.8 Number of Electrodes : The remaining determination to be made in 

the analysis is the number of electrodes required to make the 
deviations from perfect field uniformity small. To determine this, 
we simply compute the values of M in Equation (3.9) as a function of 
N, using the relations for c6 and V expressed in Equations (3.14) and 
(3.15). 

The result of such a calculation is depicted in the bar graph of 
Figure (III-4). In the graph, the RMS error (±\J"m^) is plotted for 
the cases of N = 2,4,6,8,12, and 16 electrodes. The error is 
normalized to be a percentage of V q (which is the largest potential 
appearing on the electrode boundary) and is, therefore, a plot of 
±n/V 4 versus N. 

It is noticed that the expected improvement with increasing N 
occurs for N^2, but for N = 2 there is an apparent anomally in 
this trend. 

This can be explained as follows: In the cases of N>2 (the cases 

which involve more than a single pair of electrodes) the restriction 
of requiring equal plate sizes reduces the degrees of freedom by which 
Mjj can be minimized. In other words, if the plate sizes were to 
be adjusted also, lower values of would result as compared to the 
results obtained above. In the minimization process for two 
electrodes, however, the adjustment on 06 is automatically the more 
complete variation, and so a much lower value of results. The 
equivalent restriction of equal plate size for N = 2 would be to fix 
each plate size atO(= 90° (the two plates touching) which would 
obviously give a result worse than the others, and thus, in keeping 
with the expected trend. In essence, the case of two electrodes which 
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we have calculated is not in the same class as the others. 

In fact, this case can be viewed as an N = 4 configuration, in 
which we have optimized the four plate widths as well. In so 
doing, the RMS value is halved. Even though a variation of 
individual plate widths can lead to lower values of for the 
same number of electrodes, as we have stated, it does not leave an 
electrode configuration which is amenable to field rotation techniques. 
This latter symmetry consideration is the original reason for 
placing the restriction of equal electrode widths. The important 
conclusion to be derived from this calculation is that the convergence 
of the field to a highly uniform configuration is very rapid. A 
more comprehensive way of examining the resulting field 
homogeneity is to investigate its variations inside the tank. 


3.9 Evaluation of Electric Field Uniformity Inside the Tank : In this 

section, we endeavor to determine just how uniform the field is 
inside the tank as a result of optimizing the electrode design 
parameters. We can do this by evaluating the mean square error 
integral for various 'regions of the tank using the optimum 
electrode design parameters obtained from above. We know in 
advance that the RMS deviation of the true potential from ideal 
should be much smaller near the axis of the tank than- at the tank 
walls where the potential is distorted by the discrete geometry 
of the electrodes. However, it is of interest to see how large an 
internal region possesses an RMS deviation which is below, say, 1%. 


~ 2 
To do this, we simply calculate M^ = over an area ofrfr . 

/v Q 

doing this, the expression for M^(r) iss 

oo 


Upon 


< r > - <4-2 A 1 Sin0t > • ( f > 2 + 


1 y A * ( , 

2 L_j m ( a 


.2m sin m0 6 

^ (m+1) 
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I 
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Figure III- 4 

Total RMS Error vs Number of Electrodes 
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where A. and A are the same coefficients as before, a is the tank 
radius and CO is the optimum electrode half-width. The RMS value is 
then s imp ly sj (r ) . In Figure (III-5) we have plotted the RMS 
value as a function of r for the various cases of N. The plot for 
N = 2 is the best obtainable one, in the sense of our previous 
discussions , 
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The results of this calculation show that the uniformity is 

extremely good over most of the tank volume for the cases where 

8^n£- 20. For eight electrodes, the RMS deviation from perfect 

uniformity is less than 1 % over the radial range 0£— — 0.81, 

3 

which corresponds to 67 % of the cylindrical tank volume. For 
twelve electrodes, the RMS deviation of the field from ideal 
is less than 1% over more than 80% of the cylindrical tank volume. 
For sixteen electrodes, the same homogeneity can be achieved over 
more than 987, of the tank volume. Extrapolating these results 
to 20 electrodes, definitely concludes that the RMS homogeneity for 
this case will be better than 17, over the entire tank volume. 

This case would then represent an upper limit on the number of 
electrodes required to produce sufficient field uniformity. 
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IV, THEORETICAL PERFORMANCE ANALYSIS 

4*1 Calculation of System Capacitance (Empty Tank) ; In this section 

we turn our attention to the calculation of the general N-electrode 
gage capacitance assuming empty tank conditions* 


The general procedure is to calculate Causs' integral sidttg the 

> bui. ■■■*« ■■■*■'> - 5 - ! 

line 0 = i-j; refer to Figure (IV-1). This integral gives the 

magnitude of the total charge Q per nieter residing on thi 
electrodes. The required integral is 1 , 


|i| -s f 


Cl 

/ 




dr J 


(4.1) 


where E^ is the 0 component of the interior electrifi field 
evaluated along 0 = ± By symmetry, E^ (0=~^) *? -E^ (0= so 
(1.4) becomes 


a 




jQl( ; = 2 gc I E & s(a ^) >dp 


(4.2) 


Eq can be computed from the N-electrode potential 0^(r,0) by 

! E L . 1 J&h. 

■ E e r B © s 

The form of 0^ is given by (3.7) so E^, evaluated at 0 * | f is 


(4.3) 


~a » s 




oo 

rrfS 


- " ' ‘ ' J ■/' ~ / 

_ v -.yi 2mk«- r .r, 'mjT\ 

■■•= 2 / 5 V, Cos — — ) — sin nfi Sul “r~ 

o : *s / , k N m L 

■ 3 


m = 1 k 

b in - el Iq lo e* ~c:f 

Therefore, evaluating the integral in (4.2) gives, 


5 

(4.4) 
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< V c + 2 / \ “ S ^ 


sin fh ? sin met 
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m=l 


(4.5) 
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This total charge is equivalent to the sum of the individual charges 
on each pair of plates; i.e.. 


0 = Q + 

V N v o 



to 


(4.6) 


We can now identify each of these individual changes as 


8£ V 
o o 

7 r 



mT . 

sin — sm mQ 
m 


m=l 


(4.7) 


and for each plate pair designated by k 


16 G V. 
o k 


co 


m=l 


2mk2K . . m T , 
cos — — — Sin -r- sin tna 
N 2 

m 


(4.8) 


The series appearing in (4.7) and (4.8) can be computed in closed 
form, and are 


oo 



m=l 


sin sin mQ! - , . 

2 am ^ ^ 1 -f sin a 

m 4 n l 1 - sin 0! 


(4.9) 


and 



cos 


m 


In 


1 

r,. . 

l+s in (*-jj- 4 a) 


f. . ,2k r 
l-sm(— - a) 

1 

l+s in (“^ - 0 !)] 


l-sin(— 4 - a)J 


(4.10) 


The Capacitance for each pair of plates is computed by 


.M 


2V 


(4.11) 
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Figure IV- 1 

Cross-Section of Tank Showing Path of 
Integration for the evaluation of Gauss* Integral. 


9 = 0 
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and for each k by, 

c -JiL . 

k 2V, 
k 

Thus, the N-electrode capacitance is just the parallel 
combination of C and the C ,'s or 


(4.12) 


C = C 
N o 



V 


k (4.13) 

where the upper limit on k sum is either ~or depending on 
which electrode symmetry group applies. Combining all of 
these expressions gives the general expression for the empty tank 
capacitance of an N-electrode gage. 


N 


= ir { 


In 


1-t-sin CC 
1-sin a 



In 


k=l 


l+sinC-^ + a )] [(1 - sin(^ - a )] 1 
a)] [l - sin C-l 5 ' + a)], 


l+sin(- 


/KTT 


(4.14) 


4.2 Analysis of Error Due to Fluid Configuration : In this section, 

we present an analysis of the effects of fluid Configuration on 
the system capacitance. To do this we have considered a propellant/ 
vapor configuration in which the liquid/vapor interface is § 
cylindrical surface, internally concentric with the tank wall* This 
can be made to represent two different configurations of the 
propellant when the fluid and vapor portions are interchanged (see 
Figure IV-1). This Configuration is, of course, an idealization 
but it does represent an important case because it possesses A 
symmetry which is completely immune to error reduction by orthogonal 
averaging techniques. As a result, the case is one in which the 
dielectric shielding effects cannot be reduced. 
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In the following analysis, we will derive the expression for 
the capacitance of the general N-electrode gage as a function 
of the propellant mass present in the tank. By so doing, we 
can compute the error in the measured capacitance of the 
two configurations for any amount of fuel in the tank, by 
subtracting the two results. 

We begin by solving Laplace's equation for the potential in each of 
the regions (l),and (2). The dielectric constants for these regions 
are £ ^ and 6^, respectively. The general expressions for the 
electric potentials in these regions are then given by. 


0 1 (r, 0 ) 



(4.15) 


0 2 (r, 0) 



_ (r) cos m0 
m 


(4.16) 


where, 

, » m _ -m 

<7~ (r) = a r + p r 
m m m 

(4.17) 

The task now is to determine the coefficients T , OL , and a 

m m “ m 

from the appropriate boundary conditions at the dielectric 
interface (r = b) and at the tank wall (r = a). At the dielectric 
interface we must satisfy the conditions that the tangential 
component (Eq) of the electric field is continuous across the 
boundary (at r = b) ; i.e., 


(r = b) = E q (r = b) 


(4.18) 
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and that the normal component of the displacement field (0^= €. E ) 
is continuous. Therefore, 


E„ (r = b) - -f E (r = b) 
1 1 2 


Using (4.15) and (4.16), the 0 components of the fields at the 
interface are given by 

©o 


(4.19) 


^0 1 

\< r ■ b > - - \i -Sr 


r=b 



,m- 1 . 

m yf b sin m© 
m 


m=l 


(4.20) 


and, 


E q (r = b) = 


11 ^ 0 2 
r &e 


, ,m-l „ . » (m+1). . _ 

mfa b + 6 b )sin m© 

m m 


r=b 


m=l 


(4.21) 


From these, Equation (4.18) is satisfied if 

„ -2m _ 

7- - a - p b = 0 

in m 


m m 


(4.22) 


Similarly, by computing E^ = - — from (4.15) and (4.16) 
and satisfying (4.19) gives 

'Z - ^ (a - p b' 2m ) = 0 

m 6^ m m 


(4.23) 


Now at r = a, (a, 0) must satisfy the N-electrode boundary 
condition which we calculated in a previous report. Therefore, 


tn " rn v 

a a + ^ a = <r~ (a) 

mm m 


(4.24) 
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where , 


, s 4 sin mCt 

cr^(a) = w 


m 


m 


V + 2 
o 


i 


k=l 


„ 2mkfT 
v k cos -T“ 


(4.25) 


Equations (4.22), (4.23), and (4.24) constitute the three 
simultaneous equations necessary 
£} . Solving these we obtain: 


simultaneous equations necessary for determining qr- , ot * an< ^ 


= 


2a m CT(m) 


m .* (1 +i )+b 2 - (1 .i ) 

*2 S 2 


(1 + -±) a m m) 
6 2 


a = 

m 


2m 1. , 2m 1- 

a (1 +T 4 ) + b (1 - 

z £ 2 


(1 - -i) a™ b 2m CT~ (m) 


e - ■ — ~ 

2m ^ *1. A ,2m % 

a (1 + — -) + b (1 - — ) 

6 2 € 2 


(4.26) 


(4.27) 


(4.28) 


where, as before (Equation 4.25), ^(m) = 0^(3). 


With the coefficients T , oC , and & evaluated we can now 

m m ~ m 

compute the capacitance of the gage as a function of the propellant 
mass present. This calculation is completely analogous to the one 
performed in the last report where we first determined the charge 
on the plates by evaluating Gauss' integral albng the line 
Q = ± 1^/2. Because of the two regions, the magnitude of the plate 
charge is expressed as 


Iq I = 2 


a 

e 1 j E 0 (0 = ,j |) dr + e 2 ^ E 0 (0 = dr 


(4.29) 
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where (9 = ~) and Eg (8 = are obtained from (4.20) and 
(4.21) * with 0= ~ ^ an( j b = r; 




y a - ?> - 



m=l 


m- 1 jf 

m r sm m — 
m 2 


(4.30) 


K « *?> ' 



. m-1 . n ~(m+l)v . tt 

m(Ql r +fr ') sxn m 
mm 2 


m=l 


(4.31) 

Using these expressions in (4.29) and performing the integration 
gives, after a little manipulation, 



Q I =2 \ sin m 


S 


(e r - e a )b m + y m ) + s^ b " m 

I m 2 m 2mm ' ™ 


m=l 


2 K nT 

4 

(4.32) 


From Equation (4.23) we have 


e X - e o a = - e P b 

1 . w m 2 m 2 m 


•2m 


(4.33) 


and, therefore, 


oc 


|Q | = 2 e 



, m -m. 

. sm — “r (a a -pa ) 

2 / 2 m ~ 


m 


m=l 


(4.34) 


Now, using the coefficients cc m and ft from (4.27) and (4.28) 


we have 


o® 



|q | = 2e 0 ) CT^(a) sin 


-2 ^ . 
m=l 


m 7T 
2 


m i - (b 2m ( 1 “ “)/d +T 1 ) 
. 3 € 2 e 2 


Li + (V m a - ~)/a +~) J 

2 2 (4.35) 
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where O’^(a) is given in (4.25). 


Because of the k sum appearing in J (a) we can write the total 

m 

charge Q as the sum of individual charges from each pair of 
plates; i.e., 


I - kl + Y |q i 


k (4.36) 

Using the expression for (5 m ( a ) in (4,25) we can identify each charge 
contribution from (4.35) and (4.36) and write the respective 
"geometric" capacitances* as 
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sin m'.ysin m 


zr 


m 


m= 1 
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(4.37) 
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cos ~~~ sin m?' sin mQ: 
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l_(V m (l . -!)/(l +-i) 

a e 2 e 2 


m=l 


i+(h 2m a 

a 


e € 

-r 1 )/d +“) J 

e 2 e 2 

(4.38) 


It is desirable to evaluate the trigonometric series appearing 
in these equations in closed form as we have done previously. 
However, because the coefficients in the square brackets depend on 
m we must employ a "trick" to make this evaluation possible. We 
can do this by first rewriting these coefficients in simpler form 

/- 1 (^) 2m (l-^)/d + ^) 


2m e e 

Ll + (2) (1 - ~ 1 )/(1 +- 1 ) 

a e 2 e 2 


1 - Ax 
1 + Ax 1 


m 


.m 


(4.39) 

*The geometric capacitance is the capacitance measured at the taps of 
the secondary of the transformer. However, the capacitance of 
interest is that measured on the primary side (the transformed 
capacitance) which we will compute. 
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where we have defined 


A = (1 - 


€ ! e -, 

~) / a + - 1 ) 
e 2 e 2 


(4.40) 


and 


x = (b/a) ‘ 


(4.41) 


Now we can rewrite the right side of (4.39) in a series form; 
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1 + Ax 
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= 1 + 2 
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(4.42) 


Putting (4.42) into (4.37) and (4.38) gives 
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(4.44) 
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By employing this "trick", we have broken up the original 
expressions (4.37) and (4.38) into two terms each containing 
series in m which can now be evaluated in closed fotfm. The 
evaluation of these series requires a sizeable calculation which 
we do not include here for reasons of brevity, and therefore 
simply set down the results. They are (remembering that ftt 
is restricted to odd integers only) 


©o 



m=l 


2 M u 

1 4- x + 2x sin a 
1 + x 2 ^ - 2x M sin a] 


and 


( 4 . 45 ) 
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][l+x 2 ^2x^sin(^ , '4a)] 


(4.46) 

The first terms in (4.43) and (4.44) can be obtained in closed 
form by setting x = 1 in (4.45) and (4.46). 
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and for each k. 
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(4.48) 
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The total N-electrode geometric capacitance can then be determined 
from 


N N-2 

T or — 7- 



k=l 


(4.49) 

In (4.47) and (4.48) the square of the ratio b/a is proportional 
to the fractional propellant mass present in a cylindrical tank 
when the liquid is in region (1) and proportional to one minus the 
fractional mass when the liquid is in region (2); i.e., the 
fractional mass in the first case is given by (for a tank 
cylinder of length L and radius a, and cylindrical fuel "slug" 
of length L radius b, and density p ), 


m irb Lp ,b N 2 

M = 2 ~ V 

rr a Lp (4. 50) 

In the second case we have, 


or 


2 2 

m _ (fra - 7Tb ) Lp , ,Jx 2 

M 2 . V ' 

rra Lp 


(4.51) 

Also, in evaluating the capacitance of either of the two 
different configurations of fuel we must correctly identify 

and e 2 « With these changes Equations (4.47), (4.48), and 
(4.49) can be used for both cases. 
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As we mentioned, the expressions (4.47), (4.48), and (4.49) are 
the "geometric" capacitances of the gage and do not represent the 
true measurable capacitances of the gage. This latter capacitance 
is called the transformed capacitance, and is the capacitance 
measured on the primary side of the transformer. It is 
different from the geometric capacitance because the N-electrode 
gage is viewed as an equivalent single parallel plate capacitor on 
the primary side; the geometric capacitance of each pair of plates 
is transformed differently because of the difference in the 
transformer turns ratio for each pair. If we designate the 
equivalent capacitance on the primary side as C'^, and assume that 
the transformer is ideal (loss less) and has a 1:1 voltage 
transformation at the peak voltage V q for the gage, then from 
energy considerations we can write 




~ C V 2 
2 o o 



k=l 


C k \ 


or 


(4.52) 
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C + 
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z_ . 


v k 2 

c k <r> 

o 


Now we know the voltage on each k pair of plates (from our 
previous analysis) to be 


(4.53) 


V, = V cos —~r ■ . 
k o N 


(4.54) 
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Therefore, the total measurable (transformed) capacitance of the 
gage is given by 


C' = C + 
N o 


k=l 


2 ,2kK 

cos (~f ) , 


(4.55) 


where C q and are obtained from (4.47) and (4.48). 


Incidentally, as we mentioned in the proposal, the measured 
capacitance C 1 can be enhanced (if small) by adjusting the 
turns ratio of the transformer so that the voltage on the 
primary side is less than V . 
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V. EXPERIMENTAL PROGRAM 

The experimental portion of the program was divided into two distinct 
parts. The first part consisted of a parallel work effort involved in 
the design and fabrication of the tank/ gage hardware a'nd the electronic 
circuitry. The second part involved actual experimentation with both 
the circuitry and the complete gaging system. 

Unfortunately, the experimental work, which was designed to establish 
ultimate concept feasibility, was not carried to the point where con- 
clusive evaluations of gage performance could be made. The reason for 
this is attributable to three fundamental causes. First, the experi- 
mental program was delayed about a month because of some difficulties 
which arose in the fabrication of the model tanks. Second, the program 
objectives were changed shortly after the experimental work was begun. 
Finally, the small amount of experimentation that was performed, uncov- 
ered an unexpected problem in the measurement concept, which could not 
be corrected in the allotted time. We will discuss this problem more 
fully in Section 5.5. 

5.1 Model Tank Fabrication : In designing the tank hardware, it was 

decided to keep the model relatively small so that large volumes 
of test fluids, etc,, would not have to be used in the test pro- 
gram. A dimensioned drawing of the model tank is depicted in 
Figure V-l. It was designed with two goals in mind; first, that 
in large measure, it would be a scale model of the SIVB with re- 
spect to relative dimensions (not including wall thickness); 
second, that it would be a design which would allow accessibility 
to internal surfaces during electrode installation and the test 
program. To meet the latter goal, the tank was designed with a 
flat removable top. The bottom is spun in a perfect hemispherical 
shape to simulate that portion of the SIVB geometry, and to allow 
for the investigation of electrode design in this region. The tank 
is constructed of .04 inch thick aluminum: sheet aluminum is used 

for the cylindrical portion; the end cap is spun aluminum. Two 
identical tanks were constructed. 
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Cylinder and Hemisphere Made 
From .040 TH. Aluminum 

18.50 



Fig.V-1 Missle Tank-Fuel Simulated 
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5.2 Gage Fabrication : As a result of the electrode design analysis, 

it was decided to fabricate the first gage using twelve electrodes. 
The choice is based primarily on two considerations: first, it is 

desirable to use a configuration from the orthogonal symmetry group 
N = 4, 8 12 . . . ) so that averaging could be utilized eventually; 
second, the case of N = 12 is one in which the resultant field is 
highly uniform over 807, of the tank volume, but still possesses 
some non- uniformity which can be studied. Furthermore, it was 
desirable to keep the number of electrodes as small as possible 
to ease the cabling problem and keep the electronics simpler. 

The electrode configuration was not designed so that averaging 
could be used. Thus, instead of installing patch sensors in the 
hemispherical portion of the tank, the strips were extended con- 
tinuously from one side of the tank to the other, through this 
region, each following the proper unidirectional field potential. 

The electrode fabrication required two attempts before a success- 
ful unit resulted. In the first attempt, the gage was designed 
as a separate unit so that it could be removed from tank in the 
advent of a needed replacement. This required constructing the 
electrode array on a wooden mandrel as shown in Figure V-2. The 
twelve electrodes were fabricated from copper screen of small 
mesh. The screening was held in place with tape on the mandrel 
and then coated with an epoxy resin. To prevent the epoxy from 
adhering to the wooden mandrel, the latter was impregnated with 
a sealant, waxed and coated with a thin layer of teflon. 

After the resin had cured, the outside of the electrode structure 
was machined to conform to the inside walls of the tank. At this 
point, a major problem arose. During the curing of the epoxy, 
some of the copper screening lifted away from the mandrel. This 
meant that the assembly could not be turned down enough to fit 
the aluminum tanks without machining away the screening at these 
high spots. This technique was abandoned. 
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In the second approach, it was decided to make the gage and tank 
one unit. For this case, the inside surface of the tank was coated 
with a thin layer of epoxy to provide the insulation between the 
electrodes and the tank. The electrodes were fabricated directly 
on this layer using copper foil tape. With this construction tech- 
nique, it was possible to accurately control the shape and spacing 
of electrodes. 

As part of the preparation for the test program, a test stand was 
designed and constructed which serves as a gimballed supporting 
structure for the tank. The set of gimbals allows two degrees of 
rotational freedom for the tank so that testing with liquids can 
be made at various tank orientations. The stand is depicted in 
Figure V-3. 

5.3 Electronics Design and Testing : The task to demonstrate feasibi- 

lity of the proposed capacitance measuring scheme was divided into 
two portions. First, the basic measuring scheme was demonstrated 
using an air variable capacitor to simulate the tank but without 
having the tapped transformers or the multi- segment electrodes in 
the model tank. The second portion of the task was to include the 
tapped transformers and actual tank electrodes. The reason for 
the division was to permit careful testing of the electronic mea- 
suring circuit using the precision air variable capacitor to simu- 
late the tank. The division also allowed parallel fabrication of 
the tank model to proceed. The basic measuring circuit is shown 
in Figure V-4 with both the actual and the equivalent tank circuits. 
The parameter values for the equivalent circuit are also shown in 
Figure V-4. They were chosen to be representative of the values 
to be expected from the model tank. Operating frequencies for 
initial tests are also noted. The functions of oscillator ampli- 
fier and demodulator were performed in the demonstration electronics 
by Lock- In- Amplifier is shown in Figure V-5. It can be seen that 
the reference channel tuned amplifier together with the amplitude 
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Figure V-5 

BLOCK DIAGRAM JB5 LOCK-IN AMPLIFIER 
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limited positive feedback make an oscillator. The oscillator out- 
put is fed to the external circuitry and is used as the reference 
phase for the phase sensitive detector in the Lock- In- Amplifier. 

In the signal channel preceding the demodulator, there is a tuned 
amplifier which has a characteristic similar to that of a tuned 
circuit with Q = 25. Following the demodulator there is a low- 
pass filter and DC amplifier. 

A specific experiment which has been performed will now be de- 
scribed. This experiment is the first feasibility experiment and 
used the circuit shown in Figure V-6. The first Lock- In- Amplifier , 
AR1 , operates at 100 kilo-hertz and puts out five volts peak to 
peak from its REF OUT terminals. The voltage drop across R1 , pro- 
portional to current in the primary of Tl, is applied to the SIG 
IN terminals and has a very small amount of amplitude modulation, 
synchronous with the operation of relay K1 . The output of AR1 is 
connected to the SIG IN terminals of Lock- In- Amplifier AR2 in order 
to recover the small modulation signal. The output of AR2 , then, 
should be proportional to the unknown capacitor Cl. 

The circuit was built and tested with the values as shown in Fig- 
ure V-6. These component values are those of the equivalent cir- 
cuit of Figure V-4 except that the shunt capacitors C3 and C4 are 
1,000 picofarads each instead of 3,000 picofarads. It was possible 
to obtain a full-scale output from A2 for a change in Cl of less 
than 0.6 picofarads as was planned. However, the output noise level 
was unexpectedly large (about 107 o of full scale RMS). Turning off 

the modulation completely, by disconnecting the drive to relay Kl,- 

did not reduce the noise, but disconnecting the input from AR2 al- 

most completely eliminated it. It was clear that the noise output 
of AR2 was due to noise modulation of the oscillator of AR1. A 
study of these noise sources was made. 

The actual experimental circuit shown in Figure V-6 is a simplifi- 
cation of the conceptual circuit shown in Figure V-4. In Figure V-6, 
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relay K1 has only one pole and the vacuum balance capacitor C2 
appears to the right of Transformer T1 . C2 in Figure V-6 has the 
same purpose as C2 in Figure V-4. 

Still another version of the circuit appears in Figure V-7. Here, 
the relay arm is grounded which would be desirable for an eventual 
solid-state switch. This circuit was tested and found to be much 
too dependent on stray capacitances between transformer windings. 

In the test, C2 has to be almost 90 picofarads to balance out the 
strays. While it is conceivable to use a guard shielded trans- 
former for Tl, our decision was to use the circuit of Figure V-6 
in preference to that of Figure V-7. 

The sources of noise were identified and the noise level signifi- 
cantly reduced in subsequent experiments. Principal sources of 
noise were: 

1. The High Frequency Oscillator 

2 . The High Frequency Detector 

3. Modulator Switching Transients 

The Lock- In- Amplifier was found to be noisy both in its oscillator 
and detector functions. 

The next experiment was to substitute a laboratory test oscillator 
as the high frequency oscillator. The oscillator used was a vacuum 
tube, Weir Bridge type, and was much noisier than the Lock- In- Ampli- 
fier. The laboratory test oscillator was also very microphonic. 

A switching oscillator was tried next.. This oscillator, using tran- 
sistors as switches rather than as linear amplifiers, was successful 
Its noise level was low enough so that other noise sources could be 
attacked. The circuit of the switching oscillator is shown in Fig- 
ure V-8. 

The detector portion of the Lock- In- Amplifier was also noisy; appar 
ently mostly -j noise in its output dc amplifier. A simple diode 
detector was constructed and had much less noise. 
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PROPOSED EXPERIMENTAL CIRCUIT 
FIG. V- 9b 
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The experimental circuit after the previously described improve- 
ments is shown in Figure V-9. A more detailed description of the 
circuit of Figure V-9 follows: 

The oscillator is the switching oscillator (Figure V-8). T1 is 
a 1:1: 1:1 ratio with the four windings very tightly coupled by 
quadrifilar winding. T2 and T3 are tapped onto transformers to 
give the correct electrode potentials. The taps are at 86% and 
50% of the voltage, R2 and R3 are damping resistors without which 
the tuned circuits comprised of the transformer inductance and 
shunt capacitance have too high a Q. The detector is self-explan- 
atory. The relay driver converts the sine wave at 5Hz to a square 
wave to drive the relay K1 . 

The next noise source to be attacked was that associated with 
modulation transients. The voltage across R1 in Figure V-9 drops 
to half its usual value during the "fly-time" of the relay K1 
since half the shunt capacitance is disconnected then. At the 
instant when the contact is made, there is a transient modulation 
of the envelope of the voltage across R1 . The magnitude of this 
transient depends on the phase of the high frequency oscillator 
at the moment of reconnection. The detector diode CR1 changes C4 
to nearly the peak of the envelope, thus time spreading the trans- 
ient energy. 

The next modification was to use an improved detector circuit and 
try synchronous modulation. Figure V-9b shows the new arrangement. 
The synchronous modulation is obtained by dividing down the high 
frequency oscillator to provide the low frequency modulation. The 
detector circuit was changed to one with average rather than peak 
type response to reduce time spreading of transient peaks. Also 
a blanking circuit was put in the demodulator. The blanking cir- 
cuit turns off the demodulator for a small fraction of the modu- 
lation period at each switching time, hence disregarding the tran- 
sients . 
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After incorporating these modifications, the performance of the 
system was verified using a General Radio precision variable 
capacitor. The signal-to-noise ratio was found to be greatly 
improved so that preliminary testing of the gaging system could 
begin. 

5.4 Gaging System Tests and Evaluation : During a visit by NASA early 

in December, the first gaging system test was performed in the 
form of a demonstration. 

At this time, the modifications to the electronics had been in- 
corporated, and the new tank had been completed except for some 
details of the cover; also, the fittings for filling and drain- 
ing had not been added. Accordingly, the demonstration was per- 
formed without the cover. 

Various amounts of RP-1 were poured into the tank. The results 
demonstrated a slight nonlinearity with liquid volume (approx. 

4% of full-scale). The tank was tilted in various directions 
and amounts; some sensitivity to the orientation of the surface 
was observed, but the amount of variation did not seem excessive. 

Fittings were later added to the tank to permit filling and drain- 
ing with the cover in place . 

After completing the cover and adding fittings, more data were 
taken with high resolution, and it was observed that the output 
was more nonlinear with liquid volume than had been anticipated. 

Figure V-10 shows data taken with the tank axis in two different 
orientations. There is far more separation between the two curves 
than anticipated and the curves are much more nonlinear than ex- 
pected. The reasons for these results have been analyzed, and the 
non-linearity appears to be due to three basic causes. Each cause 
produces a different non-linear characteristic in the response 
curve of Figure V-10. The three characteristics which we refer 
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to here are (1) the general separation of the two curves which is 
greater than expected, (2) the divergence of these two curves with 
increasing volume of liquid and, (3) the tapering off of the response 
at both empty and full tank. The cause for the large separation in 
curves is apparently due to a fundamental problem in the measuring 
circuit. Because it is a concept difficulty, we will discuss this 
point more fully below. 

The divergence of the curves appears to be attributable to the exis- 
tence of the center electrode which is maintained at ground poten- 
tial. Because the electrode is grounded, the fluid in the immediate 
vicinity of it does not get measured. When the tank is in the verti- 
cal position, more of this electrode is under fluid than when the 
tank is horizontal; and, therefore, less fluid is measured. This 
discrepancy increases until the tank is almost full. However, this 
difference is not large compared to the curve separation which is 
due to the concept problem. 

The tapering off of the response curves at the full and empty points 
is attributable to both the ground plane effect and the highly non- 
uniform field near the electrodes in the end caps, especially for 
the flat cover. Because the unidirectional electrode configuration 
used in this gage is not an optimum one, the existence of these 
inhomogeneities is not surprising. 

These two problems are solvable. The grounded electrode can be 
floated relative to the grounded tank wall and, thereby, be 
included in the measuring circuit by an additional set 
of leads from the transformers. The inhomogeneity problem 
can be corrected by shaping the electrodes in the dome. 
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5.5 Concept Problem : The first measurements made with a model tank. 

and breadboard electronic unit using the principles outlined in this 
report yielded two unexpected results: 

1. The change in capacitance when the model tank was partly filled 
with liquid was only about half of the expected value, 

2. The nonlinearity of capacitance change with liquid volume was 
much greater than anticipated. 

A review of the technical concepts led to the following conclusions: 

1* The concept that the use of a uniform field can be used to measure 
liquid quantity is correct; the uniform field tends to weight 
all volumes equally. 

2. The idea upon which the original electronic unit design is based 
is incorrect; the reversing of the sign of the potential on 
half of the tank boundary does not result in zero internal 
capacitance. 

The scheme outlined in the proposal and implemented in the experimental 
program has utilized a boundary potential which has alternated between 
that which caused an approximately uniform field and another with the 
same magnitude but which resulted in no current flow across a tank 
diameter. 

For the ideal case in a cylindrical tank, this would result in 
internal fields given by: 

V. (r, 0) = V - sin 0 (5.1) 

1 o a 
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and 

V (r, 9) = £ V 
l 7T o 

The purpose of this technique was to eliminate from the measurement 
the large capacitance from the electrodes to the tank walls. Because 
of the internal potential V ^ (r, 9), however, the technique also 
results in subtracting out part of the desired measurement. Since 
(r, 9) is not a uniform field, the desired equal weighting of 
each point in the region is not accomplished. 

The electric fields represented by equations (5.1) and (5.2) are 

shown in Figure V-ll. The capacitance of the field corresponding 
4 

to V„ is just — - times that of the uniform field due to V.. . 

2 n z i 

The electrodes of a uniform- field array are terminals of a multi- 
terminal capacitance network. More usefully, we may consider each 
set of electrodes which are at a single potential as a single 
terminal, thus reducing the complexity of the equivalent circuit. 

For example, consider the six electrode cylinder shown in Figure 
V-12. The potentials and appear on pairs of electrodes. 

Only four distinct potentials are used (plus ground); therefore, 
the equivalent network requires only five terminals including 
ground, as shown in Figure V-12b. Each terminal has some 
capacitance to every other terminal. 

To prevent the figure from becoming cluttered, we have labeled only 
one of these capacitances in the diagram to indicate the notation to 
be used. Network nodes are assigned letter symbols, and the capacitance 
appearing between nodes is indicated by using as subscripts the symbols 
for the two nodes involved. 


Oo 


-z 


n=i 


4m 2 - I 


(*f‘ 


cos 2ne 


(5.2) 
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As discussed earlier, the potentials applied to the electrodes 
are chosen to create, within the tank, the most uniform field 
possible with the given electrodes. The measurement required 
is then the capacitance represented by this almost-uniform field. 

At the same time, it is necessary to insure that the measurement 
is not influenced by the fields between the electrodes and 
ground. 

The capacitance which we desire to measure is readily calculated 
by considering the relationships in Figure V-12. The capacitances 
which appear directly between each terminal and ground represent 
fields which are outside of the electrode array; all others represent 
fields which are inside. 

The total electrostatic energy stored on the inside of the electrode 
array is: 

(f = 1/2 X 2Z 1/2 C jk V jk 2 j, k 4 G (5.3) 

The first 1/2 in this expression occurs because, due to the double 
sum, each capacitance is taken twice; e.g., once as and once as 

C 21* 

The capacitance we wish to measure appears across a pair of terminals 
which are driven with potential V q . The energy stored in this 
capacitance must be the same as (5.3); i.e., 

<5 - 1/2 C V 2 (5.4) 

mo 

Therefore, 

c „- (vf) < 5 - 5 > 
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Equation (5,5) states that Cm may be found by measuring the C., 

2 Jk 

individually and applying a weighting factor (V„,) to each, and 

_JiS 

adding the results; here the are the V q potentials 

which would be applied to obtain a uniform field. The actual value 
which is applied in making the measurement is, of course, immaterial. 

Equation (5,5) appears to imply that both C ^ and C 21* ^ or exam P^ e > 

are measured. Actually, of course, only one of these need be 

measured since C., = C, . . 

ik ki 

From the above discussion, it is clear that the measurement of the 
capacitance corresponding to the liquid quantity can be accomplished 
by sequential three- terminal capacitance measurements. For a system 
which has a large number of electrodes, this method may not result in 
a satisfactory system response time. 

It is equally clear that other electronic systems can be 
devised to perform the desired measurement. For example, it is 
possible to devise a switching scheme which causes the electrode 
potentials to alternate between the cosine approximation, which 
yields the uniform internal field, and a set of identical 
potentials, which yields zero internal field. By properly weighting 
the currents flowing to the various electrodes, and subtracting 
the second set of currents from the first, it is possible to obtain 
a square-wave modulation which has an amplitude proportional to the 
desired capacitance. This scheme is capable of rapid response. 
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APPENDIX A 

In this appendix a brief calculation is given for the required boundary 
potentials at the tank surfaces which give an absolutely uniform electric 
field E within the tank. Two ideal cases will be considered, namely an in- 
finitely long cylinder and a sphere. An idealized representation of the actual 
S-IVB tank (see Fig. A1 ) can be obtained by cutting the cylinder for the appro- 
priate length and fitting two hemispherical domes at each end. The bulkhead 
separating the LOX tank from the LH 2 tank is also assumed to be hemispherical, 
for convenience. This idealization results in a spherical LOX tank. The tank 
problem is solved by equating the potentials and electric fields along the 
junctions . 

We consider the cylinder first. Since we assume the cylinder is of in- 
finite length the problem is simply two dimensional; see Fig. A2 . The appro- 
priate coordinates for the problem are the polar coordinates (r, 9). It is 
easy to see from the figure that if the electric field E is to be everywhere 
constant and uniform within the circular boundary r = R and in the direction 
indicated, then the line AA' is a ground potential (V = 0), by symmetry. Thus 
the product of the field E and the perpendicular distance from AA' to the point 
(r, 0) must be equal to the potential difference V(r, 0) between these points. 

Therefore , 

V(r, 0) = E r cos 0. (A-l) 

Similarly the potential at the boundary r = R must be 

V(R,0) = E R cos 0. (A- 2) 
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We can arbitratily define ER = V Q , so that the potential anywhere in- 
side r = R becomes 

V(r, 0) = V_ £ cos 9, (A- 3) 

° R 

and the boundary potential is dependent only on 9 through the relation 

V(9) = V 0 cos 0 (A- 4) 

Therefore, on the cylindrical portion of the ideal tank we must apply the 
potential (A-4) to produce the uniform field E = V Q /R everywhere inside. 

In the case of the LOX tank, which is assumed to be spherical for conven- 
ience, the problem is three dimensional. Again the cross-section AA' (see 
Fig. A3 ) is at ground potential by symmetry and thus the same field E times 
the perpendicular distance between AA' and the point (r, 0, 0) must be equal 
to V(r, 0, 0). Thus inside the tank, 

V(r, 0, 0) = V 0 — sin 0 cos 0, (A-5) 

R 

and at the wall, r = R we have 

V (R, 0, 0) = V Q sin 0 cos 9. (A-6) 

In these relations we have defined V Q = ER, as before, so that the values of 
E and V at the junction of the hemisphere and cylinder are matched. 

The potential distribution (A-6) completely solves the cases of the 
simplified geometry of the LOX tank and the domes of the hydrogen tank, 
while equation (A-4) solves the problem of the cylindrical LH 2 tank wall. 
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APPENDIX B 

The Mean Square Error Integral Function 


The minimization of the Mean Square Error Integral Function is a standard 
mathematical technique for representing a known function in terms of a linear 
series expansion of orthogonal functions. To take a simple example, which is 
analogous to the problem in the text, let us assume that we want to represent 
an arbitrary function f (x) , as well as possible, over an interval a^x£f b. 
Furthermore, we want to represent it by a series sum of a complete set of 
orthonormal functions 0 k (x) defined over the same interval. Such a set 
might be sines and cosines or Bessel functions, etc. This is immaterial. 
Mathematically, we desire: 


f(x) = aL A k 0 k (x); 


over a 




6 b 


(B-l) 


where the A k are constant coefficients. To represent f(x) accurately by the 
right hand side, we must adjust the magnitude of each of the A k appropriately. 
In general, at each point x there will be some discrepancy, or error, e (x) 
in this equality which we can represent as: 


f(x) = ^ A k 0 k (x) + € (x) 


or 


e(x) = f(x) - X/ A k 0 k (x) 


(B-2) 


We can reinterpret f(x) as an average value of ^ A k 0 k , and the problem 

then becomes one of adjusting the A k so that e(x) is a very small deviation 

at each point. This amounts to saying that we want ^ A k^k to as c -*- ose 

to its average value f(x) as possible. The procedure for minimizing the 

deviations e (x) is just one step removed from just computing the statistical 

RMS deviation of a set of data from its average; i.e., if one has a large set 

of N data points X- , with an average X = — Xj x^ then the average squared 

1 N i 

deviation is, by standard statistical procedure: 


(Ax) 2 



(B- 3) 
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The Mean Square Error Integral Function defined in the text is just 
this statement written for a continuous set of data (functions) instead of 
a discrete set. The correspondence is: 

.2 


z 


J 


x 


x- 


> 


M 

b , 

. . .dx 

a 

f(x) 


A k 0 k (x) 

k 

b 


i.e., different "data points" 
for different A k> 


N > (" ^ dx 

J a 

Written out for the continuous case' using B-2, the analog becomes. 


M - « 


[e(x)J 2 dx f b jf(x) ”^A k 0 k (x) 

v- a 


‘dx 


dx 


dx 


j 


a 


(B-4) 


Now, the fitting of S A k 0 k (x) to f(x) is equivalent to minimizing M by 

k 

adjusting the "data points" so that they coincide with the average; this 
is equivalent to a dishonest experimentalist adjusting every x^ in (B-3) 
so that each is equal to "x (the number he wants) and, thus, the data has 
no experimental errors. In our case, we just adjust the A k so that M is 
a minimum and we can do this by a calculus of variations technique; i.e. ; 
by computing: 


J5M 

3A k 


= 0 


(B-5) 


In the present example, if carried out using B-4, this would result in an 
expression for the desired A k , which is: 

A k = _ 1 _ P f(x)0 k (x)dx 
b-a J a 

Notice that if the 0 k were sines or cosines with the interval endpoints 
b =iT, a = 0 then this expression would be the familiar Fourier coefficients 
relation; i.e.. 


r TT 


A k = 


rr 


f(x)cos x dx 


'->0 


In terms of the uniform field problem, we are using this procedure (which is more 
involved computationally, but the same in principle) to find the coefficients V k 
and tX. which allows 0N(<^, r >©) to represent 0 O (r,Q) as closely as possible. 


-B2- 



